Converting Zitterbewegung Oscillation to Directed Motion 



O 



Qi Zhang, 1 Jiangbin Gong, 2 > 3 <EI and C.H. Oh 1 ' 4 '!! 

1 Centre of Quantum Technologies and Department of Physics, 
National University of Singapore, 117543, Singapore 
2 Department of Physics and Centre for Computational Science and Engineering, 
National University of Singapore, 117542, Singapore 
J NUS Graduate School for Integrative Sciences and Engineering, Singapore 117597, Republic of Singapore 
^Institute of Advanced Studies, Nanyang Technological University, Singapore 639798, Republic of Singapore 

(Dated: January 21, 2013) 

Zitterbewegung oscillation (ZB), namely, the jittering center-of-mass motion predicted by free- 
space Dirac (or Dirac-like) equations, has been studied in several different contexts. It is shown 
here that ZB can be converted to directed center-of-mass motion by a modulation of the Dirac-like 
equation, if the modulation is on resonance with the ZB frequency. Tailored modulation may also 
stop, re-launch or even reverse the directed motion of a wavepacket with negligible distortion. The 
predictions may be examined by current ZB experiments using trapped-ion systems. 
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I. INTRODUCTION 

Zitterbewegung motion (ZB) originally refers to the 
trembling motion of a free relativistic particle (e.g., elec- 
tron) as described by the Dirac equation [1]. ZB in this 
context has an extremely large frequency and a very 
small amplitude, so it is hardly observable. Neverthe- 
less, the physics of ZB becomes highly relevant to a num- 
ber of interesting situations where the quantum dynam- 
ics is governed by artificial Dirac or Dirac-like equations. 
Examples include band electrons in graphene [2|, cav- 
ity electrodynamics Q , ultracold atoms in designed laser 
fields I^HSli macroscopic sonic crystals and photonic su- 
perlattices @, as well as trapped-ion systems [l(| [HI- 
For a review of recent fruitful studies of ZB physics using 
electrons in semiconductors, see Ref. [12]. 

Many quantum dynamical phenomena induced by a 
control field have been explored using the Schrodinger 
equation. Much less is known in the case of driven Dirac 
or Dirac-like equations. This work is concerned with ZB 
in a simple version of a driven Dirac-like equation. Pre- 
viously, we showed that ZB in the cold-atom context as 
a quantum coherence phenomenon can be coupled with 
some mechanical oscillations, resulting in the control of 
the amplitude, the frequency, and the damping of the ZB 
associated with a wavepacket Q. Going one step further, 
it should be of interest to study what happens to ZB if 
one term in a Dirac-like equation is modulated with a 
frequency on-resonance with the ZB frequency. Indeed, 
in both classical mechanics and quantum mechanics, it 
is always expected that an oscillation phenomenon can 
behave remarkably different when it is subject to an on- 
resonance driving field. A latest example is Bloch oscilla- 
tion subject to on-resonance driving, where Bloch oscil- 
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lation may be converted to 
or even to directed motion 

In this work we start from a Dirac-like equation that is 
already experimentally realized in the context of trapped- 
ion systems [ll[ • We then introduce a modulation to the 
Dirac-like equation, whose frequency is on resonance with 
the ZB frequency. Due to the on-resonance driving, it is 
shown, both analytically and numerically, that ZB as- 
sociated with a wavepacket can be converted to directed 
wavepacket motion. By tailoring the modulation, we may 
also stop, re-launch or even reverse the directed motion 
of a wavepacket with negligible distortion. These results 
indicate that ZB is not a mysterious oscillation, and it 
provides an extra useful control knob for manipulating 
quantum motion. 



II. THEORETICAL ANALYSIS 

Let us consider a simple version of a two-component- 
spinor Dirac-like equation (taking h = 1 throughout), 

ij t \iP)=H\iP) = ~(Ca y +Dp x a x M), (1) 

where p x is the momentum in the x direction; C and 
D are two c-number coefficients, o x ,<J y are Pauli ma- 
trices in the x and y directions. We first assume C is 
time-independent, with C = Co > 0. Time-dependence 
of C will be discussed much later. We also note that 
the Hamiltonian H here depends on p x , but not on the 
momentum along the y direction. Such kind of Hamil- 
tonian may be synthesized by tripod-scheme cold atoms 
with a strong constraint potential in the transverse direc- 
tion Q, tripod-scheme cold atoms in a two-dimensional 
laser setup with mirror oscillations [8j, and trapped-ion 
systems Qij, HH ■ Because the a x p x term is the only spin- 
orbit coupling term, the ZB oscillation associated with a 
wavepacket does not suffer from a fast damping in its 
oscillation amplitude @, Q . 
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We now examine the motion of an initial state (t = 0) 
as a product state of a one-dimensional wavepacket in x 
and an internal two-component spinor state, 



cos 2 (C t/2) sin 2 (C t/2) 



<a#(0)) = G(x) 



^Pox 



(2) 



where G(x) is a broad Gaussian in real space centered 
at xo = 0; po (assumed to be zero below) is the central 
momentum of the wavepacket along x. To better under- 
stand ZB, the dynamics can be analyzed in the momen- 
tum space. Hence we carry out the Fourier transforma- 
tion of the wave function to get a narrow wavepacket in 
momentum representation (the required narrowness will 
be explained later), 



(3) 



where g(p x ) is also a Gaussian as the Fourier transfor- 
mation of G{x). 

With the initial state specified in Eq. the in- 
ternal state will evolve because it effectively experi- 
ences two "magnetic fields": one along y of strength 
Co and the other along x with strength Dp x for the 
component p x [see Eq. (1)]. The total effective mag- 
netic field strength is yCg + Px-D 2 and the direction 
of the total magnetic field is characterized by an an- 
gle arctan(Z?p x /Co). To obtain analytical results, we 
keep effects up to the first order of Dp x /Co, thus mak- 
ing approximations aictaii(p x D/Ca) ~ p x D/Co and 
\/Cq +p x D 2 « Co- Physically, our approximation is to 
assume that for different p x components, their effective 
Zeeman splitting is assumed to be almost the same, but 
with the internal state precessing around slightly differ- 
ent directions characterized by arctan(Dp.j;/Co). In par- 
ticular, the analytical solution to the Dirac-like equation 
([T]) can be directly obtained (for C = Co), i.e., 

m)) = cos Qc t) 9(Pz) ( I ) e- ix ° p » 

+ sinQc t)^ x )^)e-^ + W^. (4) 

Here an extra independent phase factor e^^ ^ 00 ^^ in 
the second term arises from the p x - dependence of the 
spinor precession axis. The meaning of our early assump- 
tion of a narrow wavepacket in the momentum space is 
also clear: the main part of the wavepacket should ensure 
that \p x D\ <C Co- Considering effects due to (p x D/Co) 2 
or high orders will capture some insignificant deformation 
of the wavepacket and the damping of ZB, but as demon- 
strated by our full numerical results below, these high- 
order effects will not change the essence of the physics 
under discussion and can be neglected here. 

Interestingly, the analytical solution in Eq. (fj| can be 
interpreted as a superposition of two sub-wavepackets 
centered at Xq = and at x = Xq + D/Cq = D/Cq- The 
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FIG. 1: Schematic plot of wavepacket motion satisfying a 
Dirac-like equation. In case (a), a Dirac-like equation is 
not modulated. Cyclic population transfer between two sub- 
wavepackets gives rise to ZB oscillation. In case (b), a Dirac- 
like equation is periodically modulated, with the modulation 
frequency on-resonance with the ZB frequency. Due to the on- 
resonance modulation, sub-wavepackets successively reemerge 
in different positions, yielding directed motion. 



occupation probabilities of the two sub-wavepackets are 
given by cos 2 (Co£/2) and sin 2 (Cot/2). As time evolves, 
these two occupation probabilities oscillate, thus giving 
rise to a time-dependence of the average position de- 
noted (x). The jittering motion is between x = and 
x = D/Co, with the angular frequency Co- Specifically, 
(x) is given by 

(x) = cos 2 (C t/2)xo + sm 2 (Cot/2)(x + D/Co) 

= 2^[l-cos(Cb*)], (5) 

which is nothing but a ZB phenomenon, with the ZB 
frequency Co, the ZB period Tzb = 27r/Co, and the 
"equilibrium point" of the oscillation at x = D/(2Co). 
In particular, when t = ±tt/Cq, the sub-wavepacket lo- 
cated at x = disappears and the other one located at 
x = D/Cq reaches its maximal population. Figure 1(a) 
schematically depicts this ZB mechanism. 

Given this interpretation of ZB in terms of the pop- 
ulation cycling between two sub-wavepackets, we next 
investigate what if the sign of the a y term is changed at 
t = tt/Cq = Tzb/2, with the modified Hamiltonian given 
by 



H = lj(-C (Ty + Dp x a x ). 
At that instant, the state of the system is given by 



iPx\ip(t = tt/Cq)} = g(p x ) 



A 



1 



(6) 



(7) 



Adopting the same approximations as introduced above, 
we analytically obtain 

m)) = cos[C (; - Tr/C )/2]g(p x ) ( J \ 

+ sin[C (t - 7T/C )/2]g(p x ) ( J ) e -* ( W^(8) 
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Remarkably, Eq. © represents a superposition of a sub- 
wavepacket located at x — D/Cq and a new one at 
x = 2D/Cq. The new equilibrium position of ZB is 
hence shifted to x = 3D/(2Cq). Loosely speaking, this 
indicates that if a sign change of the Ca y term occurs 
at the right moment, then the /^-dependent extra phase 
factor [see Eqs. (4) and (8)] in the time-evolving state 
makes an opposite contribution and consequently a sub- 
wavepacket at a new location emerges. From this new 
solution, it is also observed that at t = 2ir/C = T ZB , 
the occupation probability of the sub-wavepacket cen- 
tered at x — D/Cq will totally vanish and that of the 
other centered at x = 2D /Co will become unity. The 
average position of the system moves from x — D/Cq at 
t = Tzb/2 to x = 2D /C at t = T ZB . 

Repeating this strategy, i.e., abruptly changing the 
sign of C after every interval of Tzb/2, a directed trans- 
port emerges from a modulated ZB. That is, if we let 



(a) 2T / T = even 

stop ZB 



C 



mod(; 



(9) 
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, 1), then ZB physics no longer in- 



where /; = 
duces oscillations. Instead, ZB is forced to become di- 
rected motion, as schematically illustrated in Fig. 1(b). 
Certainly, modulating the Ca y term in a slightly differ- 
ent manner may also result in directed transport in the 
— x direction. We will not repeat the similar theoretical 
analysis. In addition, as shown below, tailored modula- 
tion of C can lead to more complicated control over the 
wavepacket dynamics. 



III. NUMERICAL SIMULATION 

In this section, we verify our analytical results by nu- 
merical simulations without any approximation. Two 
representative modulation schemes are considered. In 
scheme (i), we introduce the modulation of C as de- 
scribed by Eq. © first and then stop the modulation for 
Tstop = NTzb [N integer), followed by the same mod- 
ulation again [see Fig. 2(a)]. In scheme (ii), everything 
is the same except that 2T st0 p/IzB is an odd integer [see 
Fig. 2(b)]. 

For scheme (i), Fig. 2(c) shows that the wavepacket 
first moves in the x direction, then undergoes ZB once 
the modulation is stopped, and finally the directed mo- 
tion continues after the modulation is switched on again. 
The wavepacket profile in real space at different times is 
plotted in Fig. 2(e). For scheme (ii), Fig. 2(d) shows 
that the wavepacket first moves in the same way as in 
the first scheme. However, once the modulation is re- 
launched, the wavepacket transport is reversed. So this 
behavior may be regarded as a type of "super-ZB" oscil- 
lations (analogous to super-Bloch oscillations [Hi). As 
also shown in Figs. 2(e) and 2(f), in both schemes the 
distortion of the wavepacket profile for the shown time 
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FIG. 2: (color online) Numerical results for a modulated 
Dirac-like equation. Panels (a) and (b) illustrate two schemes 
of modulation of the sign of system parameter C. When mod- 
ulation is on, the sign of C is reversed after each time interval 
of Tzb/2 [see Eq. {SJ]. The modulation off-period r st0 p is an 
integer multiple of Tzb in (a) and a half-integer multiple of 
Tzb in (b). Panels (c) and (d) depict the time-dependence of 
the expectation value of (x), for modulation schemes in (a) 
and (b), respectively. It is seen that in (c), the wavepacket 
continues its directed motion after the modulation is switched 
on again; but in (d), the wavepacket eventually reverses its 
directed motion. Panels (e) and (f ) show the wavepacket pro- 
file, i.e., probability density distribution in x, at three differ- 
ent times, for modulation schemes in (a) and (b), respectively. 
It is seen that wavepacket distortion is hardly visible. Note 
that in (f), the wavepacket at t = 125.6 is on top of that at 
t = 0. All plotted quantities and system parameters are in 
dimensionless units, with Co = D = 1. 



scale is invisible to our naked eyes. The numerical results 
here confirm our theoretical analysis, demonstrating that 
ZB can be indeed converted to directed wavepacket mo- 
tion in a well-controlled fashion. 

In real experiments, it should be more straightfor- 
ward to introduce a smooth modulation to system pa- 
rameters. Consider then a sinusoidal modulation of C, 
C = Cosin(o;i). This case is somewhat complicated be- 
cause the magnitude of C, and hence the ZB frequency, 
is also time-dependent. To resolve this technical issue, 
we define a time-averaged effective ZB angular frequency 
over one half period of modulation, i.e., 



J ZB — 



7r/u 



i=0 



2C 

C sm(ujt)dt= — -. (10) 

7T 



In the light of our previous interpretation of ZB oscilla- 
tion, it is expected that after the time interval 7r/w|g, 
one of the two sub-wavepackets has vanished and the 
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FIG. 3: (color online) Same as in Fig. 2, except for a sinusoidal 
modulation of C, i.e., C = Co sin(ajt), with uj = 1, Co = vr/2, 
and D = 1. The pattern shown in panels (c)-(f) here are 
much similar to those seen in Fig. 2(c)-2(f), but the distance 
traveled by the wavepackets is different. 

other has an occupation probability of unity. In accord 
with our theoretical insights above, the sign of C must be 
changing at this moment (from CP to + or vice versa) 
in order to convert ZB oscillation to directed motion. 
This suggests the following modified resonance condition: 
uj = = 2£o. It should be also noted that for a si- 
nusoidal modulation, the magnitude of C can be very 
small during certain time windows and hence our early 
treatment in terms of the first-order expansion of D/C 
can be somewhat problematic. Nevertheless, as seen in 
Fig. 3, the numerical results confirm that oj — cj|g is the 



required resonance condition under a sinusoidal modula- 
tion so that ZB oscillation is again converted to directed 
motion. This also implies that the controlled dynamics 
is mainly contributed by the time segments during which 
C is appreciably nonzero. As seen in Fig. 3, for either 
even or odd 2T stop /TzB (Tzb is now determined by 
the pattern of the results presented in Fig. 3(c)-(f) is es- 
sentially the same as those seen in Fig. 2(c)-(f), with 
again little wavepacket distortion. Note also that here 
the wavepacket transport distance per modulation cycle 
is less than that achieved by a mere discontinuous sign 
modulation of C. 



IV. DISCUSSION 

The Dirac-like equation in Eq. ([T]) can be realized in 
cold-atom systems or in trapped-ion systems. Accord- 
ing to the concrete experimental realization in Ref. fTlj j . 
the Ca v term can be realized through two on-resonance 
laser fields, with the system parameter C connected with 
the phase and the strength of the laser fields. Modula- 
tion of C can hence be realized by modulating the laser 
phase, via a laser phase modulator or even an oscillating 
mirror Q. The ZB frequency in Ref. [ll[ can be sev- 
eral tens of kHz, and the required phase modulation fre- 
quency should match this ZB frequency. Again using the 
experimental parameters in Ref. , it can be estimated 
that for a ZB frequency of 40 kHz, a wavepacket can be 
transported by about 10 -8 m per modulation cycle. 

In summary, using a simple Dirac-like equation that 
is already experimentally realized, we show that ZB can 
be converted to directed motion via a periodic driving 
of the Dirac-like equation, if the driving is on resonance 
with the ZB oscillation. We hope that our results can 
motivate further studies of quantum control in driven 
Dirac-like equations. 
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